Back-action effects in an all-optical model of dynamical Casimir emission 



O 
(N 

m 



Oh' 
I 



> 

p 
o 



% 



I. Carusotto/'G S. De Liberate,^ D. Gerace,^ and C. Ciuti^ 

^ INO-CNR BEC Center and Dipariimento di Fisica, Umversita di Trento, 1-38123 Povo, Italy 

^Laboratoire Materiaux et Phenomenes Quanttques, 

Universite Pans Diderot- Pans 7 and CNRS, UMR 7162, 75205 Pans Cedex 13, France 

'^ CNISM-UdR Pavia and Dipartimento di Fisica "A. Volta", Umversita di Pavia, 27100 Pavia, Hah 

(Dated: July 5, 2011) 

We report a theoretical study of the optical properties of a three-level emitter embedded in an op- 
tical cavity including the non-rotating wave terms of the light-matter interaction Hamiltonian. Rabi 
oscillations induced by a continuous wave drive laser are responsible for a periodic time-modulation 
of the effective cavity resonance, which results in a significant dynamical Casimir emission. A clear 
signature of the back-action effect of the dynamical Casimir emission onto the drive laser is visible 
as a sizable suppression of its absorption. 
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One of the most intriguing predictions of modern quan- 
tum field theory is the possibility of converting the zero- 
point fluctuations of a quantum field into real particles 
when the boundary conditions of the field are varied in 
time at a fast enough rate. The most celebrated exam- 
ple of such an effect is the emission of electromagnetic 
radiation from a moving but electrically neutral metallic 
mirror, the so called dynamical Casimir effect (DCE) [l|. 

So far, experimental observation of the DCE has been 
hindered by the extremely weak intensity of the emission 
in realistic configurations |2| , but theoretical understand- 
ing of the underlying physics can be considered as almost 
satisfactory. The next theoretical challenge consists in 
combining the quantum field theory with a microscopic 
description of the mirror motion, including the mechani- 
cal back-action effect of the quantum field onto the mir- 
rors: pioneering theoretical works have anticipated the 
appearance of a very weak friction force to compensate 
the radiated energy [3| , but to the best of our knowledge 
no experimental investigation has been reported yet. 

New perspectives in the study of the dynamical 
Casimir effect were opened by the realization that a 
sizable DCE can be obtained with a suitable time- 
modulation of the effective optical length of the cav- 
ity: many proposals to very the optical length have been 
put forward, from photo-generation of a highly reflecting 
plasma in the mirrors [J, l5| , to time- modulations of the 
refractive index of the cavity material [6| , to coupling the 
cavity mode to an emitter with time-dependent proper- 
ties [jjW,- First experimental evidences of a DCE have 
recently appeared using a superconducting circuit termi- 
nated by a SQUID [1^; as theoretically proposed in \U\ . 
the effective electric length of the circuit is modulated by 
driving the SQUID with a time-dependent magnetic flux. 

In this Letter we consider a model of a cavity whose 
effective optical length is modulated in time in an all- 
optical way by means of a drive laser beam that in- 
duces Rabi oscillations in a three-level emitter in a lad- 
der configuration strongly coupled to an optical cavity. 



as sketched in FiglTJa). Such a configuration has been 
recently realized in the experiment of jl3| using a com- 
bination of inter-band and inter-subband transitions in 
a semiconductor device and a theoretical investigation of 
its optical properties at the level of the rotating wave 
approximation (RWA) has appeared in [13|. Of course. 



the same model is applicable to other physical realiza- 
tions, such as atoms in microwave cavities [ij], Joseph- 



15|, or quantum 



son qubits in superconducting circuits 
dots in photonic crystal cavities |la |. 

Here we go beyond the RWA and perform a complete 
quantum optical calculation of the system dynamics in- 
cluding all the relevant non-RWA terms in the light- 
matter interaction Hamiltonian as well as in the dissi- 
pation superoperators. In particular, we anticipate the 
appearance of an appreciable quantum emission by dy- 
namical Casimir effect. More importantly, the optical 
nature of the cavity modulation allows to read out the 
back-action effect of the quantum emission from the ab- 
sorption that the drive laser experiences upon interacting 
with the emitter. This effect is the optical analog of the 
friction force experienced by the moving mirrors in stan- 
dard DCE configurations J4I, but its actual measurement 
is expected to take great advantage of the extreme pre- 
cision of optical techniques over mechanical ones. 

The internal dynamics of the emitter-cavity system is 
described by the model Hamiltonian 
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where fkOg^ej are the energies of the g,e, f atomic lev- 
els, respectively. The <? -H- e transition is optically driven 
by a coherent laser of frequency lol and (real and posi- 
tive) Rabi frequency ^eg- The standard RWA has been 
performed on this transition under the assumption that 
its frequency Wg ~ ^g is much higher than all other fre- 
quency scales. A single cavity mode of frequency Wcav 




FIG. 1: (Color online) Left (a) panel: sketch of the ladder 
level configuration in the emitter coupled to the cavity and of 
the optical fields driving the transitions. Right panels: total 
emission intensities 7cav and J/e as a function of the Rabi 
frequency Q^g of the drive laser (solid lines) . Red dashed line 
in (b): same curve in the vanishing emitter-cavity coupling 
case f^cav = 0. Within the RWA both /, 
zero. The chosen parameters oje 
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correspond to a strong (but not ultra-strong) coupling regime. 



and destruction (creation) operator a (a^) is considered, 
which is strongly coupled to the e — >■ / transition with a 
(real and positive) vacuum Rabi frequency flcnv In order 
to correctly describe the DCE, all terms of this coupling 
have to be taken into account, including the anti-RWA 
ones where a cavity photon is emitted while the atom 
climbs from the e to the / state and viceversa. In the 
following, we shall restrict our attention to the resonant 
case with uj^ = uJe — ujg and ojcav = i^/ — Wg. As re- 
quired by the strong light-matter coupling regime, the 
light-matter coupling J7cav is assumed to be much larger 
than all decay rates. 

Thanks to the RWA assumption on the g ^ e transi- 
tion, we can apply the unitary rotation operator R(t) 



transition is assumed to be forbidden), as well as via 
direct light emission from the cavity through the non- 
perfectly reflecting mirrors. Such dissipation processes 
are included in the model at the level of the master equa- 
tion for the density matrix p, 

^=-^[H,p]+ C,g [p] + Cf, [p] + £eav [p] ■ (3) 

Thanks to the large value of We — ujg, the standard RWA 
approximation can be performed on the e — )■ 5 spon- 
taneous emission terms, which leads [17| to a super- 
operator which has the usual Lindblad form, /Ceg[/o] = 
h^eg[2(T-gP(T+g-(T+ga-gp-p(T+g(T-g] lu both thc statlouary 
and rotating frames, with (t+ = |e)(g| and a^ = (''^eo)^- 
Nonetheless, it is worth reminding that energy conser- 
vation shows peculiar features when seen in the rotating 
frame of the dressed Hamiltonian ^■. as the emission 
spectrum on the e — )■ g transition is shifted downwards 
by ll)l, the sidebands on the red side of Wi = uje~ujg have 
a negative energy in the rotating frame; on the dressed 
level scheme of Figl2l[a,d), these transitions appear in fact 
as climbing up the energy ladder. 

Consistency of the theoretical model requires includ- 
ing the anti-RWA terms for the spontaneous emission 
on both the / o e transition and in the cavity emis- 
sion. This requires some additional care as a number of 
counter- intuitive optical processes may occur where, e.g., 
the atom goes from the e to the / state while emitting a 
photon by spontaneous emission. For a sufficiently weak 
coupling to the baths [18], the dissipation superoperators 
can be written in the temporally local form 

^j [p] = r, {u,ps, + Sjpu] - SjUjp - pu]Sj}, (4) 

where the system-bath interaction operators for j — 
{/e,cav} have the form Sfe — |e)(/| + |/)(e| and S'cav — 
a^ + a, respectively. Energy conservation in the decay 
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Hamiltonian has a time-independent form 
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The structure of the resulting eigenstates of the cavity- 
emitter system optically dressed by the drive laser is 
shown in Figl21[a,d) for different values of ^eg ■ The labels 
of the eigenstates refer to the dominating component in 
the weak flcav hmit, in which the energy of the \gn) ±\en) 
and \fn) eigenstates are nwcav T ^eg and {n + l)ajcav re- 
spectively. In the shorthand |jn), n and j — {g,e,f} 
respectively indicate the number of cavity photons and 
the state of the emitter. 

In addition to the cw driving laser, the system is cou- 
pled to the external world via spontaneous emission pro- 
cesses on the e —>■ g and f ^r^ e transitions (the f —?' g 



whose definition involves the Fourier transform of the 
frequency-dependent density of states of the baths, 
Vj{t) = (27r)~^ /^dwe"'""^ Wj(w). As a result, the de- 
cay from the |in) to the |fin) eigenstate of the dressed 
system occurs at a rate Tj Kfinj^jlin)!^ ^(win — L^Hn)- In 
the present paper, all baths are assumed to be zero tem- 
perature ones, so that energy can only dissipated from 
the system into the bath, which imposes v(uj) = for 
a; < 0. For the numerical calculations shown in the fig- 
ures, we have considered model forms of Vj{uj) such that 
Vj{uj) = for either w < or cj > cjmax, with a UV 
cut-off Wmax choscn to be much higher than all other en- 
ergy scales of the problem. In the intermediate region 



< U) < UJr, 



(w) is taken to be flat and equal to 1 



exception made for a smoothening of the edges. We have 



checked that ah physical conclusions do not critically de- 
pend on the specific choice of the UV cut-off Wmax and 
of the details of the smoothening 22 1. 

The expectation value of all one-time operators of the 
system in the stationary state are immediately obtained 
from the steady-state solution pss of the master equation 
([3]), in particular the full photon number distribution in 
the cavity. However, as it was pointed out in [19|, this is 
not enough to determine the actual light emission from 
the system: to distinguish real radiation from the virtual 
photons that are present even in the ground state of the 
cavity, one has in fact to evaluate the full spectral distri- 
bution of the emitted radiation [l7| , that is the Fourier 
transform of the two-time correlation function |23l] 

/oo 
dTe^-^{S,{t + T)S,{t)) (6) 
-oo 

for j — {cav, /e}, from which the total emission inten- 
sities Ij are obtained after integration over w. Plots of 
Icavje as a function of fleg are shown in Fig[Hb,c). 
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FIG. 2: (Color online) Left (a,d) panels: sketch of the dressed 
levels of the optically driven emitter as predicted by the 
Hamiltonian ((2|). The labels indicate the dominant contri- 
bution of each eigenstate in the weak atom-cavity coupling 
ficav limit. Right panels: spectra of the spontaneous emis- 
sion on the / -fT- e transition (b,e) and of the cavity emission 
(c,f). Upper and lower group of panels (a-c) and (d-f) refer 
to the Q,eg/u]ca.v = 0.7, 2 cases, respectively. Same system 
parameters as in Fig[l] 

Both emissions have a negligible value [2J| for low Rabi 
frequencies J7eg < Wcav/2. The first threshold occurs at 



structure in /cav for larger f^eg's, in particular two peaks 
at rigg — cjcav and 2a;cav When the emitter is not cou- 
pled to the cavity (il^av — 0, red dashed line in FiglTJb)), 
the lower threshold and the two peaks disappear, while 
the upper threshold is almost unaffected. On the other 
hand, the emission completely vanishes if the anti-RWA 
terms are not included. 

A physical explanation of these numerical observations 
is obtained by looking at the level schemes of the opti- 
cally dressed system that are shown in Figl2l^a,d). For 
Q.eg < a;cav/2, the dynamics of the system is limited to 
the subspace spanned by the two lowest energy eigen- 
states \g, 0) ± |e, 0) at energies ^Hileg- For fleg > ^cav/2, 
the energy of the \gn) — \en) dressed states starts exceed- 
ing the energy of the \g{n + 1)) -I- |e(n + 1)) state, which 
activates a novel family of decay processes indicated by 
the blue solid arrows in FiglJJa). These decays occur via 
spontaneous emission on the e — > / transition, the matrix 
element being provided by the weak admixture (propor- 
tional to ficav/^eg) of the neighboring \fn) state into the 
final \g{n + 1)) + \e{n + 1)) state by the emitter-cavity 
coupling ficav As a result, population is transferred to 
the upper manifolds and significant emission intensities 
Icavje appear as the dressed system goes back to the 
lowest manifold. The second, stronger threshold that is 
visible in Fig[TJb,c) at fleg = i^cav is due to an additional 
e — ^ / spontaneous decay channel that opens up as soon 
as the \gn) + \en) state exceeds in energy the \fn) state. 

This interpretation of the thresholds is confirmed 
by the plots of the frequency-resolved emission spectra 
G'cav(w) and Gfe{oj) shown in Figl2jb,c). Note that, in 
contrast to the previous discussion of the e — >■ g emis- 
sion, the frequency w in Gcav,/e(w) corresponds to the 
actual physical frequency of the emitted light. The cen- 
tral peak at w ~ Wcav corresponds to transitions between 
dressed states that only differ by the number n of cav- 
ity photons. The lateral peaks originate from transitions 
between different dressed states, e.g. at cj ~ 2fleg — Wcav 
{\gn) — \en) -^ \g{n + 1)) -|- |e(n -I- 1)), solid blue arrows) 
and u! ~ Wcav ±^eg ilfn) — )■ \gn) ± |en), dashed and dot- 
ted red arrows). Correspondingly to each arrow in the 
diagram (a), the vertical lines in the spectra (b,c) indi- 
cate the expected position of the peak: the agreement 
with the numerical spectra is excellent. 

Similar reasonings can be used to explain the peaks 
that are apparent in Fig[TI^b,c) around fleg = Wcav, Swcav^ 
the stronger one at Q^g — 2a;cav originates from the reso- 
nant mixing of the \gO) — |eO) and |/1) states by the anti- 
RWA terms of the emitter-cavity coupling (see the hori- 
zontal arrow in Fig[2jd)). The peaks at Wcav ± f^cav/V2 
in the frequency-resolved spectra of Figl2je,f) indeed cor- 
respond to transitions to and from the new eigenstates in 
the form of linear combinations of the \gO) — |eO) and |/1) 
states. The interpretation of the weaker and narrower 



,/2. A second, more pronounced threshold is peak at Vie 



in Figlljb,c) is similar: the mixing 



now occurs between the \gO) — |eO) and the \g2) + |e2) 
states, the weaker matrix element being due to the smah 
non-resonant admixture of |/1) to the \g2) + |e2) state. 
The peak around fleg — 2tJcav has a simple physical in- 
terpretation in terms of the dynamical Casimir emission: 
under the effect of the driving laser, the atom performs 
Rabi oscillations on the g -^ e transition, so that the 
atom-cavity coupling is periodically switched on and off 
at frequency 257eg tlSj. As usual, amplification of the 
zero-point fluctuations is most effective when the modu- 
lation frequency is close to an even multiple of the cav- 
ity frequency [3|- The radiated energy is compensated 
by those e ^ g spontaneous emission processes that 
move the system upwards in the dressed level scheme 
of FiglJJd), i.e. the spontaneous emission of photons at 
a physical frequency slightly smaller than the one of the 
absorbed photons from the drive laser at uj^. 




FIG. 3: (Color online) Solid line: Rate of photon absorption 
from the drive laser as a function of its Rabi frequency fleg. 
Same system parameters as in Figd] Red dashed line: same 
curve in the vanishing emitter-cavity coupling case ficav ~ 0. 
Blue dotted line: same curve for a two level emitter. 

We finally turn to the back-action of the emitter onto 
the drive beam. The simplest effect to consider is the ab- 
sorption of the drive laser light by the emitter performing 
the g — >■ e transition. Using standard results of quantum 
optics [17[ , the average rate of photon absorption is re- 
lated to the imaginary part of the expectation value of 
the emitter polarization by Reg — 2fleglm{Tr[cl pss]}- 
A plot of this quantity as a function of fl^g is shown in 
Fig El For low fteg, the dependence is the standard one 
of saturated absorption by a two-level atom (blue dot- 
ted line). The downward jumps at n^g — '^cav/2 and 



(more visibly) at fie 



are closely related to the 



thresholds previously observed in Icav and Ife- 

Clear signatures of the back-action effect by the dy- 
namical Casimir emission are visible as absorption dips 
around ileg = uJcav and fleg = 2wcav If the frequency of 
Rabi oscillations is close to the resonance peak at which 
dynamical Casimir emission is maximum, the emitter in- 
deed experiences a sizable friction force that tends to op- 
pose the Rabi oscillations and effectively reduces the rate 
of energy absorption from the drive [25| . This interpre- 
tation is validated by the observation that both negative 



peaks completely disappear for ficav = (red dashed line 
in Figl3]) . Remarkably, the amplitude of both peaks is a 
non-negligible fraction of the total absorption rate. 

In conclusion, we have theoretically studied the dy- 
namics of an optically driven three-level emitter in a lad- 
der configuration embedded in a single-mode microcav- 
ity. This has required taking into account those novel 
optical processes that stem from the anti-rotating-wave 
terms of the light-matter coupling Hamiltonian. Rabi 
oscillations under the effect of a continuous wave driv- 
ing laser result in a sizeable dynamical Casimir emission 
from the cavity. A clear signature of the backaction of the 
dynamical Casimir emission onto the source of the time- 
modulation is visible as a marked dip in the absorption 
from the drive beam when the dynamical Casimir emis- 
sion is the strongest. Generalization of this approach to 
more complex geometries such as optical black holes [20[ 
is expected to shine light on fundamental problems of 
contemporary physics such as black hole evaporation un- 
der the effect of the Hawking radiation 21 1 . 
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